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Abstract. Recent rotation measure observations of a dozen or so galaxy 
clusters have revealed a surprisingly large amount of magnetic fields, 
whose estimated energy and flux are, on average, ~ 10^^ ergs and ~ 10^^ 
G cm^, respectively. These quantities are so much larger than any coher- 
ent sums of individual galaxies within the cluster that an efficient galactic 
dynamo is required. We associate these fields with single AGNs within 
the cluster and therefore with all galaxies during their AGN phase. Only 
the central, massive black hole (BH) has the necessary binding energy, 
~ 10®^ ergs. Only the accretion disk during the BH formation has the 
winding number, ~ 10^^ turns, necessary to make the gain and magnetic 
flux. We present a model of the BH accretion disk dynamo that might 
create these magnetic fields, where the helicity of the a — dynamo is 
driven by star-disk collisions. The back reaction of the saturated dy- 
namo forms a force-free field helix that carries the energy and flux of the 
dynamo and redistributes them within the clusters. 



1. Introduction 

The problem of understanding the origin of the large scale galactic magnetic 
fields has been with us for over forty years. There have been many papers 
and reviews on the galactic and extragalactic magnetic fields (see Moffatt 1978; 
Parker 1979; Krause and Radler 1980; Ruzmaikin et al. 1988; Wielebinski & 
Krause 1993; Beck et al. 1996; Zweibel & Heiles 1997; Kulsrud 1999), and 
observational reviews (see Miley 1980; Bridle & Perley 1984; Kronberg 1994), 
including the observations themselves (e.g. Perley et al. 1984; Taylor et al. 
1990; Taylor & Perley 1993; Eilek et al. 1984). 

Recent rapid progress in observational work on galaxy clusters has revealed 
a surprising result. Intracluster medium (ICM) appears to be definitely magne- 
tized, and in many cases, perhaps are highly magnetized as convincingly argued 
by Eilek et al. (2000). Figure 1 presents one such example in Hydra A Cluster 
as shown by the rotation measure (Rm) map made by Taylor & Perley (1993). 
We will show in this article that the implied magnetic energy and flux estimated 
from extensive Rm maps of a dozen or so galaxy clusters are so exceedingly large 
that the conventional galactic dynamo models may prove to be inadequate. We 
argue that a new source of energy and a different form of the galactic dynamo 
are required. 

As the rotation measure observations of galaxy clusters are relatively new 
and some of them are (yet) unpublished by the observation teams, we will first 
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explain some of the observation results in detail, then discuss their physical 
implications at length. In the second half of the article, we will propose a new 
paradigm related to AGN accretion disks and describe some of our recent efforts 
in understanding a sequence of physical processes revolving around the origin of 
cluster magnetic fields. 



2. Galactic and Extragalactic Magnetic Fields 

Faraday rotation measures, Rm, are shown to be consistent with six other phys- 
ical interpretations of magnetic fields in ours and nearby galaxies (star light 
polarization, interstellar Zeeman splitting, synchrotron emission, synchrotron 
polarization, and inferred by x-ray emission and cosmic ray isotropy and pres- 
sure) (see Kronberg 1994 for a review), thus establishing as a reliable mea- 
sure of galactic and extra galactic magnetic fields. Because of the existence of 
many self-illuminating as well as background sources, usually AGN, and the in- 
creasing sensitivity of radio detection, Rm has become the recognized measure 
of extragalactic magnetic fields (Kronberg 1994; Taylor et al. 1994; Ge &: Owen 
1994; Krause & Beck 1998). 

2.1. Magnetic Flux and Energy in Galaxy Clusters 

Recently, high quality Rm maps of self-illuminating sources of galaxy clusters 
where the distances are known have become available (for example, Taylor &; 
Perley 1993; Eilek et al. 2000). An important quantity that has received less 
discussion in these papers is the magnitude of the magnetic flux and energy. 

Figure 1 shows the Rm, map of the region illuminated by Hydra A in the 
cluster (courtesy of Taylor k. Perley 1993). The largest single region of highest 
field in this map has approximately the following properties: the size L ~ 50 kpc 
and B ~ 33)U G, derived on the basis that the field is patchy and is tangled on a 
4 kpc scale. This leads to a startling estimates of flux, F BL^ ~ 5 x JO^^/x G 
kpc?, and energy, W = {B"^ /8tt)L^ ~ 4 X 10^^ ergs, assuming that the tangled 
field is only confined to the 50 kpc region. If this is extended to the whole cluster 
which is ~ 500 kpc, then the implied flux and energy are correspondingly larger 
by a factor of 100 and 10'\ respectively. A similar conclusion can be reached 
when a larger sample of Rm of galaxy clusters are analyzed using the data 
presented in Eilek et al. (2000). In Table 1, we have reproduced part of table 
given in Eilek et al. (2000) and added two columns where the approximate flux 
and energy are calculated assuming that the fields are partially tangled or in 
loops. 

Furthermore, the estimated values of fluxes and energies are most likely to 
be the minimum of the actual magnetic fields existing in the galaxy clusters. 
Faraday rotation depends upon the component of the field strength along the 
line of sight, B\\, the distance along the line of sight, Zg, and electron density, 
Ug. Estimates of can be made from the x-ray emission measurements of the 
clusters with a typical accuracy of ~ 20%, and it varies by factors of 2 to 4 
over the region of the source, but otherwise is nearly uniform, and clumping is 
small (Taylor et al. 1994). If the field is folded in any fashion so that regions of 
oppositely directed field are in the line of sight, then the observed Rm will be 
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Table 1. A list of cluster core parameters and their estimated mag- 
netic fluxes and energies. Mean magnetic field is taken to be ~ 
\/3 X Data in the first three columns are taken from Eilek et 
al. (2000). 



Source 


Size 

(kpc) 


(^11 ) 


(1058 ergs) 


(lO^^G cm2) 


A 400 


100 


2.9 


3 


5 


A 1795 


7 


18 


0.03 


0.17 


A 2052 


8 


17 


0.05 


0.17 


A 2029 


10 


1 


0.0003 


0.02 


A 2199 


30 


15 


2 


0.3 


A 2634 


140 


1.9 


30 


7 


A 4059 


10 


69 


0.15 


1 


CygA 


70 


15 


25 


1.5 


Hydra A 


50 


33 


40 


15 


Virgo A 


3 


35 


0.01 


0.05 



smaller than that if the same field lines were straightened out into one direction. 
In other words Rm is a minimum measure of -By. 

To put the above numbers in perspective, for a typical galaxy like ours, 
e.g., with 1 kpc thickness, 3 kpc Romberg radius, and a field of ~ 3^ G, the 
magnetic flux and energy are roughly 10'^^ G cm^ and 4 x lO"''^ ergs, respectively. 
One observes that the flux and energy given in Table 1 range from close to the 
Hydra A limit to no more than 10^ times that of a typical galaxy. 

The magnitude of the implied fluxes and energies are so large, xlO^ and 
xlO^ respectively, compared to these quantities within standard galaxies that 
their origin requires a new source of energy and a different form of the dynamo 
than previous galactic models. These minimum energies are sometimes even 
larger than the baryonic binding energy of galaxies (~ 2 x 10^^ ergs). The 
extremely large fluxes also seem out of reach via amplification by ordinary galaxy 
rotations in a Hubble time. 

Next, we discuss the difficulty with using turbulence to create these nearly 
uniform, highly correlated and coherent regions of Rm as seen in Figure 1. We 
then discuss the still greater difficulty of creating the total magnetic energy of 
the cluster based upon a turbulence dynamo model. 



2.2. Turbulent versus Coherent Fields 

It has been suggested by a number of people (Eilek 1999; DeYoung 1980; Ruz- 
maikin et al. 1989; Goldman &; Rephaeli 1991; Jaffe 1980) that the entire cluster 
is uniformly turbulent due to Rayleigh Taylor instabilities during matter in-fall 
into the cluster, and that this turbulence drives the cluster dynamo creating the 
flelds. The problem with this interpretation is the total magnetic energy, the 
magnitude of the turbulence, the strength of the fields, the apparent correla- 
tion of R„i maps with single AGN structures, and finally the limited number of 
rotations of the cluster in a Hubble age. 
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Because of the small rotation rate of the typical cluster, ~ 100 km s , 
the available rotation energy is small, ^ 10~^ of the cluster binding energy, 
which has a thermal velocity of ~ 10^ km s~^. So applying the turbulent model 
to Hydra A implies a magnetic energy 10^ greater than the rotational energy. 
Therefore the dynamo must be of the type where fields are generated on the 
small scale, yet as Taylor et al. (1994) point out, the fields of Hydra A and 
A1795 reverse on the different sides of the core, requiring coherence on scales of 
100 kpc. Since this reversal is correlated with the structure of the source, and 
since the energy generated at the small scale is small compared to the turbulence 
input and the turbulent input should be small compared to the binding energy 
(DeYoung 1992), we believe that all these factors point to random, localized 
sources of magnetic energy of size > lO^'' ergs. This is probably too demanding 
for turbulence. 

Furthermore, it will be difficult to produce large scale coherent Rm regions, 
which have been observed in Hydra A (Figure 1, northern region) and several 
other galaxy clusters (Eilek et al. 2000) . This is because in a turbulent plasma, 
the emission, the R„i, and the degree of polarization should all be statistically 
symmetric. Despite the unlikelyhood of all these factors conspiring to create 
both a pattern and a nearly uniform R^, one observes in many Rm maps of 
AGN, mostly in clusters, a distinct match in the Rm pattern with the jet like 
pattern of emission. Particularly the sign of the average Rm in several cases 
reverses across a symmetry plane through he core of the AGN (Taylor &i Perley 
1993). The size of the regions of uniform Rm correlates strikingly with the size of 
the jet as a function of distance from the nucleus. We interpret this correlation 
as due to the source of the field being the AGN jet as opposed to a turbulent 
a — Q dynamo in the cluster as a whole. 

2.3. Average Field Structure 

Using serendipitous polarized background sources, therefore random lines of 
sight through random clusters, Clark & Kronberg's (1999) have made ~ 80 
Rm measurements. Their observations have produced a boundary of the typical 
cluster in Rm such that the average field was 3^ G out to a radius of Rduster — 
300 kpc. The magnetic fiux and energy, 10^ /Lt G kpc? and lO^'^ ergs, are then 
similar to the largest structure already discussed in Hydra A. If each galaxy of 
a typical chistcr with ~ 50 large galaxies contributes a high field region during 
its AGN phase, then the probability of intersecting such a region of area that 
is ~ 1% of the cluster is roughly c± 5 x 10~^ so that in 100 lines of sight, the 
probability of intersecting a Hydra-like region of an AGN in a cluster is 50%. 
This is not inconsistent with the variability they observed. Finally we note that 
the large degree of polarization observed ~ 50% in these sources indicates that 
the rotation source and emission source cannot be in the same location (Burn 
1966; Taylor 1991), otherwise polarized emission from various depths in the 
source would undergo different degrees of rotation and hence emerge depolarized. 
Therefore in any model the Faraday screen and the emission source must be 
related and even congruent in order that the screen and hence Rm be correlated 
with the core of the AGN. 
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2.4. Black Hole Accretion Disk as the Engine 

Purely based on the energetics, the accretion disk around supermassive black 
holes in AGNs offers an attractive site for the production of magnetic fields. 
The accessible binding energy of the black hole is ~ ergs 
and the winding number of the disk forming the BH of nearly every galaxy 
is ~ 5 X 10^° at WRg, where Rg is the BH horizon (r^ lAU). Using the 
canonical numbers thought to apply to AGN disks, the BH dynamo flux can be 

-pRHdyn - ^BHdynTri^BHdyn-^^' - "^^^"^ ^ Cm^, where We have used SeHdyn - 10^ 

G at Lagn ~ 10^^ ergs s^^ and i?BHdyn — lOi?^ ~ 10^^ cm. Both the flux 
and energy from this simple analysis are ~ 10 times the maximum observed 
values. No other source of energy is likely to be sufficient by many orders of 
magnitude. Therefore it is much more reasonable to assume that every AGN, 
both within and external to clusters, produces the magnetic energy and flux 
that we observe in this extreme case from the binding energy released in the 
accretion disk forming the central BH. This implies that every galaxy contains 
a BH where ~ 90% — 95% of the accessible binding energy is transformed into 
magnetic energy during its AGN phase by an accretion disk dynamo. On the 
average this flux and energy is distributed throughout the universe as force-free 
fields and only a small fraction 5% — 10% of the magnetic energy is dissipated in 
the form of the AGN spectra, thus explaining the problem of the missing AGN 
luminosity (Richstone 1998; Krolik 1999). In this picture a larger fraction of 
the magnetic energy is dissipated where the brightest AGNs are seen in galaxy 
clusters, because only in the clusters is a sufficient gas density retained by the 
gravity of the cluster such that this density confines the field increasing the 
fraction of the magnetic energy that is dissipated. For most galaxies external to 
dense clusters a small fraction of this magnetic energy is dissipated as the AGN 
radiation, a small fraction remains in the galaxy, and the bulk of the magnetic 
energy and flux is distributed in the walls and voids of the universe. 



3. Astrophysical Requirements and Progress with a Model 

The sequence of phenomena that can explain this astonishing extragalactic mag- 
netic flux and energy must start with an accretion disk forming a massive central 
galactic BH. This in turn presumes an answer to an equally enigmatic question, 
namely the formation of these massive galactic BHs themselves (Begelman et al. 
1989; Rees 1999). By focusing on the transport of angular momentum we believe 
that the flat rotation curve mass distribution can be explained as a plausible re- 
sult of any non-linear collapse of an initial gaseous baryonic density fluctuation 
by hierarchical tidal torquing (Newman &; Wasserman 1999). The BH forms 
from this mass distribution when the Rossby vortex torque mechanism super- 
sedes tidal torquing and an accretion disk forms. All this mass then collapses to 
a BH. The fiat rotation curve, M on R, results in S cx R"^. When this thickness 
reaches E ~ 100 to 1000 g cm~'^, heat is confined for several revolutions, and 
the Rossby vortex instability initiates at M^ig^ ^ 10'' to 1O*M0. Finally the 
dynamo produced fields then supersede the previous torque mechanisms. 
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3.1. The Rossby Vortex Torque Mechanism 

We have predicted and demonstrated analytically and numerically how a new 
instability in Keplerian flow, the Rossby vortex instability, can grow (Lovelace 
et al. 1999; Li et al. 2000a; Li et al. 2000b). The production of vortices is 
shown in Figure 2. This instability produces torque and thus transports an- 
gular momentum within an accretion disk by purely hydrodynamically means 
via the interaction of large, two-dimensional, co-rotating Rossby vortices. The 
enhanced transport of angular momentum by co-rotating vortices is recognized 
in rotational atmospheric flows (Staley & Gall 1979) and in laboratory mea- 
surements of the Ranque-Hilsch tube (Hilsch 1946; Frohlingsdorf & Unger 1999, 
Colgate & Buchler 1999). 

3.2. The Dynamo, Star-Disk Collisions, and Helicity 

A coherent dynamo can form in a Keplerian accretion disk because of the large 
azimuthal velocity shear, provided that there exists a robust source of non- 
axisymmetric helicity. Classically turbulence has been invoked to explain this 
helicity using the mean field dynamo theory, but we know of no way to create this 
degree of turbulence, with vertical motions, hydrodynamically in an accretion 
disk, because hydrodynamic turbulence alone is damped in an accretion disk 
(Balbus & Hawley 1998). The magnetic instability of Balbus &: Hawley will lead 
to turbulence, but the magnitude of the turbulence is orders of magnitude too 
small compared to the Keplerian stress. Instead wc have identified a new, robust 
source of helicity driven by star-disk collisions by a small mass fraction ~ 10~^ — 
10"^ of pre galaxy-formation stars. The Keplerian shear and a star-disk collision 
with the twist producing helicity is shown in Fig. 3. We have demonstrated by 
laboratory flow visualization experiments of how plumes, driven in a rotating 
frame, counter rotate relative to the frame (Beckley & Colgate 1998; Beckley et 
al. 2000) and thus produce a robust and coherent helicity where flux is always 
added in the same direction and where the driving force is large compared to 
the Keplerian stress in the disk. 

Wc have simulated the positive, exponential gain of both the quadrupole 
and dipole poloidal field of such a dynamo with a vector potential code in 3-D, 
cylindrical coordinates, where the velocity field simulates both the Keplerian 
rotation and star collision-produced plumes. We have observed a growth rate 
of ~ 10% per revolution, two plumes per two revolutions, Rpiume = ^/•^Rdisk, 
and with a magnetic Reynolds number, Rgy qb = 100 (Pariev, Colgate Sz Finn 
2000). 

3.3. The Saturation of the Dynamo and the Formation of the HeUx 

With positive gain and large winding number, the dynamo will saturate re- 
gardless of how small the seed field is. Since the helicity does not depend on 
turbulence, it will not be subject to turbulent a-quenching at the small scale 
(Vainstein &: Cattanio 1992; Vainshtein Sz Rosner 1991). Furthermore since the 
stars maintain virial velocity, their velocity is supersonic relative to the disk and 
the resulting shock stress is large. At the back reaction limit, the field grows until 
the torque of the field affects the Keplerian motion, and the accessible BH bind- 
ing energy is converted into magnetic energy. The progressive loss of this flux is 
a force-free helical, Poynting magnetic flux, which we identify as the collimated 
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AGN jets. We have investigated the field topology of these twisted helical flux 

surfaces by integrating the Grad-Shafranov equations for a force-free axisym- 
metric field with a Keplerian distribution of winding number (Li et al. 2000c) 
as shown in Fig. 4. Since the field decreases as B^gu^ oc l/R, the pressure, 
at large radius as the helix extends to Mpcs, becomes of the order of the IGM 
and and the outer boundary of the helix is self collimating (Lynden-Bell 1996). 
The energy carried by this helix, at a mean radius near the BH, Rdyn — ^ORbh 
is the accessible energy of accretion or Mbhc^/W = {B^^i-^/87r){100'!TR'^jj) or 
Bhelix - 10^ G, = 5RheiixBhelix = 5 X 10^^ amperes, and Vpotential = 10^° 
volts and / x Vpotential = watts = 10^^ ergs s"^. General relativity inside 
the innermost stable orbit will add additional energy (Blandford & Znajek 1977; 
Livio et al. 1999). 

3.4. J|| Reconnection and Acceleration 

The distribution of this flux in the universe occurs by partial tearing mode 
reconnection producing the minimum energy Taylor state (Taylor 1986). The 
total flux is conserved, but a fraction of the energy is dissipated in the tearing 
mode, J|| reconnection. The resulting acceleration of the current carriers 
produces the emission that we associate with the AGN. 
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